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Abstract 

We consider several gauge invariant higher dimensional operators that couple gravity, 
gauge fields and scalar or fermionic fields and thus break conformal invariance. In particular, 
we consider terms that break conformal invariance by the photon coupling to heavy and 
light fermions. While the coupling to heavy fermions typically do not induce significant 
magnetic fields, the coupling to light fermions may produce observable magnetic fields when 
there are a few hundred light fermions. Next we consider Planck scale modifications of the 
kinetic gauge terms of the form f{(f))F^jjF^^ and h{il))F^jjF^^'^ , where / and h are functions 
of scalar and pseudoscalar fields 0, ip^ and F^^, F^^" and the gauge field strength and its 
dual, respectively. For a suitable choice of / sufficiently strong magnetic fields may be 
produced in infiation to be potentially observable today. The pseudoscalar coupling may 
lead to birefringence in infiation, but no observable magnetic field amplification. Finally, we 
show that the photon coupling to metric perturbations produces by far too weak fields to be 
of cosmological interest. 
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1 Introduction 



Inflation H is at this moment the only paradigm that offers explanation of structure forma- 
tion consistent with current microwave background 0, |^ and large structure observations, 
as well as the observed homogeneity and isotropy of the Universe on large scales and the 
flatness, horizon and monopole problems. In this paper we consider several mechanisms in 
which conformal invariance of the gauge fields in inflation is broken and discuss under what 
conditions inflation can produce magnetic fields on cosmological scales. These fields may of- 
fer an explanation for the seeds of the micro-gauss magnetic fields observed in spiral arms of 
galaxies today 0] , but also for the intergalactic magnetic fields that might influence the cos- 
mological microwave background radiation (CMBR) and early structure formation. At the 
moment there is only an upper bound on cosmological magnetic fields of about 10~^ Gauss 
based on the current CMBR measurements In this paper we argue that one can obtain 
observable cosmological magnetic fields from inflation. Hence, magnetic fields should be 
included into the analyses of the upcoming high precision CMBR measurements, in partic- 
ular in the measurements that include polarization. To support this point of view we note 
that there is now evidence for intergalactic micro-Gauss magnetic field [Q, and for magnetic 
fields of strength 10^^ Gauss correlated on 1 Mpc from statistical correlations between 
the positions of high-energy cosmic rays and BL Lacertae sources (quasars with jets). 

It is well know that quantum loop corrections to general relativity induce corrections to 
Einstein's gravity [§] in curved space-time backgrounds. Based on this observation Turner 
and Widrow |^ have argued that breakdown of conformal invariance in inflation may be 
responsible for large scale magnetic fields that are potentially observable. In particular 
they considered (i) the gravitational coupling TZA^^A^, TZ^^A^A^, that induce a photon 
mass at the price of breaking gauge invariance, (ii) gauge invariant dimension six terms 
W^fTf^^Fp^/m^ n^^gP^F^pF^Jm^, TZg^"" gP"" F^^Fp^/m'^ , and (iii) couplings of the photon 
field to other fields (scalars and axions). Here TZ^^p^ , TZ^'^, TZ denote the curvature tensors and 
scalar, respectively, F^^, = d^Ay—dvA^ denotes the field strength, and m is a mass parameter. 
The result of their analysis is that Type (i) models may lead to strong magnetic fields today, 
because the photon couples directly to the curvature, which is large in inflation, TZ ~ H^, 
and sufficiently small today so that it does not conflict with the current observational bounds 
on the photon mass. The analyses of Type (ii) and (iii) models were inconclusive. 
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In Ref. |Tn] Dolgov has shown that the trace anomaly due to the couphng of gauge fields 



to (massless) fermions may lead to gauge field amplification in infiation. Ratra [|TT] has 
considered the effect of the infiaton coupling to gauge fields of the form e~^'^^^^'' g^^g"'^ F^j^Fp^^ 
and showed that one can obtain large scale magnetic fields in extended infiation. This idea 
was then reconsidered in the context of string cosmology by Gasper ini, Giovannini and 
Veneziano [|l^ and by Lemoine and Lemoine |]l3l, where the role of the infiaton is taken 



by the dilaton field. Garretson, Field and Carroll |]14[ have considered coupling of the 
electromagnetic field to a pseudoscalar boson; their conclusion was that it is difficult to 
obtain large scale magnetic fields in this mechanism. Mazzitelli and Spedalieri [|15] have 



reconsidered Type (ii) models of Ref. in the light of the Schwinger-DeWitt expansion |16 



for a charged scalar field and argued that one cannot obtain large magnetic fields from 



these interactions. Calzetta, Kandus and Mazzitelli [0 have considered the charged scalar 
current in scalar electrodynamics as a source for large scale magnetic fields; their results have 
been contested by Giovannini and Shaposhnikov |T^. Magnetic field production from the 
dynamics of extra dimensions has been considered in |19[. Finally, Maroto |^ has recently 
suggested that coupling of gauge fields to metric perturbations may lead to significant gauge 
field amplification. 

In Ref. I^TI we have shown that the backreaction of superhorizon scalar fields breaks 
conformal invariance of gauge fields in infiation, leading to a gauge field generation mecha- 



nism. Alternatively, the photon may acquire a mass in infiation by a Higgs mechanism [22 



when a charged scalar field gets an expectation value driven by a negative coupling term, 
~ — to the infiaton field s. When applied to electromagnetism, this then results 
in magnetic field generation with the spectrum Bi oc where i is the correlation length. 



The resulting field is sufficiently strong ||21|, |23|, ^ to seed the galactic dynamo mechanism 
in fiat universes with dark energy component. The galactic dynamo may be the mechanism 
that explains the micro-Gauss magnetic fields observed in many galaxies today [Q. An alter- 
native is the Biermann battery which may be operative on megaparsec scales at galaxy 
formation ||26|| . 

Magnetic fields may also be generated at cosmological phase transitions P? 



as a conse- 



quence of hypermagnetic field amplification in presence of the right handed electrons through 
the Abelian triangle anomaly [EB, EH] , and in presence of a pseudoscalar field condensate [BD 
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These causal mechanisms create magnetic fields which are correlated typically on too small 
scales to be of relevance for galactic and super-galactic magnetic fields today For a 

review of other models for cosmological magnetic field generation see . 



In this paper we focus on mechanisms for magnetic field generation from terms induced 
by gauge field couplings to other quantum fields in presence of gravity. We restrict ourselves 
to terms that break conformal invariance in gauge invariant manner. The paper is organized 
as follows. In section |^ we consider the evolution of gauge fields is presence of heavy fermions 
in a curved space-time background. We make use of the local Drummond-Hathrell action |P 
and find that magnetic field amplification in infiation is typically very small, unless the 
number of massive fermions Np is very large, that is Np ^ 10^. Next, in section 13 we 



consider the anomalous gauge field coupling to light fermions and find a potentially 
observable amplification of gauge fields in infiation, provided the number of light fermions 
is of the order a few hundred. In section ^ we then consider how modification of the gauge 
kinetic terms —^F^^F'^'' and —\F^yF^^ through the coupling to scalar and pseudoscalar 
fields may affect the dynamics of gauge fields in infiation. We find that the scalar coupling 



may result in significant amplification of gauge fields |TT|. On the other hand, the coupling 
to pseudoscalar fields leads to an interesting example of birefringence, but results in no 
observable amplification of magnetic fields. We then in section ^ show that gauge field 
coupling to scalar cosmological perturbations |^D[ can lead only to a very minute amplification 
of gauge fields in infiation. In section || we briefiy discuss preheating, and section ^ contains 
conclusions. 

2 Photon coupling to gravity with heavy fermions 
2.1 Drummond-Hathrell action 

We now reconsider the effect of coupling of gauge fields to gravity, corresponding to Type (ii) 
models of Ref. [Q), and argue that one cannot obtain strong magnetic fields from infiation. 
Drummond and Hathrell have shown that the scatterings of photons on gravitons in 
presence of massive fermions with a mass M (shown in diagrams on figure |l|) induce a 
correction to the Einstein-Hilbert action which, when computed in the expansion. 
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leads to the following effective action 



K 1 

a^^a^'^F F + £i 



S = J d x^— ^ 

Ci = -^[b'JZg^Pg''''F,,Fp^ + cn^''gP''F,,F,^ + dn^''P''F^,F, 



pcr\ 



where G = (SnMp)^^ denotes Newton's constant, /?e = ae^F/lSOyr, = e^/47r, where e is 
the photon coupling constant at the scale relevant in inflation, Np the number of fermionic 
species with mass M, b = —5, c = 26 and d = —2, TZ denotes the curvature scalar, TZ'^'^ and 
J^^lupa curvature tensors, F^;^ = S^Aj^ — the gauge field strength, g = det[5'^,y] the 
determinant of the metric g^i^, and M the fermion mass. In case when the gauge field is 
non-Abelian with an SU(N) symmetry, Np should be replaced by Np — llN/2 [|T^]. Coupling 
to non-Abelian gauge fields thus reduces Pe and, as we shall show, damps the gauge field in 
inflation. This is contrary to fermionic loops, which induce amplification of gauge field. The 
action (0) is computed in the expansion, and hence should be relevant for the photon 

field dynamics only when the curvature scale, which is characterized by the Hubble parameter 
H, is small in comparison to the mass scale, \TZ\ ~ <^ M^. Moreover, the physical 
photon momentum fcphys should satisfy A^phys ^ l"^!- When \7l\ -C fcphys ^ the term 
— {Ql3e/M'^)g'^PDp{g^'^Ffj_r,)g"^D^Ffjy, induced by the one-loop vacuum polarization diagram 
in figure |^, may become relevant. For large physical momenta, k'^T^y^ ^ M^, the 
expansion in (|1]) breaks down, and the fermions should be treated as massless. Dolgov |]T0[ 
has shown that in this case the fermion loop in figure ^ contributes to the stress-energy trace 
anomaly, which we consider in section 3. The problem of magnetic field growth in inflation 
in terms of the Schwinger-DeWitt expansion |TB[ for scalar fields has been studied in ||15|| . 



We shall now study the dynamics of gauge fields given by the Drummond-Hathrell ac- 
tion (|l]) in conformally flat space-times. Recall that inflation, radiation and matter era are 
all conformally flat space-times, whose metric can be parametrized as g^j^i, = a?{T)iri^ij, where 
Vpv = diag(l, —1, —1, —1) denotes the Minkowski metric, and the scale parameter a = a(r) 
is a function of conformal time r. 

The equation of motion for the time-like field Aq is nondynamical, and can be easily 
inferred from equations (|l]) and (p6| ) : 



VMo -drV-A = 0. (2) 
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Figure 1: The scattering diagrams that contribute to the Drummond-Hathrell action. The gravi- 
tons {G^j^y) are represented by the dashed blue hnes, the photons {A^) by the wavy red lines and 
the fermions (/) by the solid green hnes. 




Figure 2: The one- loop diagram contributing to the higher order in derivatives dimension six 
operator of the Drummond-Hathreh action. 



which can then be used to remove Aq from the dynamical Aj-equation. The resulting equation 
is purely transverse, and identical to one obtained in Coulomb gauge, Aq = Q = V ■ A. In 
Appendix A we show that the dynamical photon field equation can be recast as (c/. Eq. (|101|) ) 



dr 



1 + 2/?e^ (26 + 33m;) j drAj 



2/?e^(10 + 3^)] V^Af 



0, 



(3) 



where w = p/ p, (3e = agNp /180tt. In case when there are fermions with different masses, 
/3e/M^ should be replaced by where Pi = aeNpi/lSOTr. For example, when the 

number of charged massive fermions is of the order Np ~ 10^, which is to be expected in a 
typical grand-unified theory, and taking ~ 1/30, we then have (3e ~ 1/15. Since typically 
H/M <^ 1, the terms that break conformal invariance are small. We shall now make this 
statement more quantitative. 



2.2 Photon field in infiation 



Note first that in de Sitter inflation {w = p/p 



-1) Eq. @ becomes conformally invariant 



M2 



T 



0, 



(4) 



which is in agreement with Ref. j^. Hence, in de Sitter inflation there is no photon field 
amplification. For power law inflation however conformal invariance is broken. Consider 
now the homogeneous mode, for which equation is solved by 

drA[ = Co 1 + 2/5e— (26 + 33w) , (5) 



where Cq = Cq{x) characterises the initial electric field in inflation. For w < —26/33 the 
field drA^ grows in inflation. When 2/3e(-ff^M^)(26 + 33w) ^ — 1 a large amplification may 
result. In this case adiabaticity is broken for super horizon modes, leading to amplification. 
To make this more quantitative, we now quantize the field as follows 



E 

T=l,2 



Ak-x 



—ik-x 



(6) 



where A^^ {i = 1, 2) denote the gauge field mode functions, e| (T = 1, 2) are the transverse 
polarization vectors, and and ^ the photon field annihilation and creation operators. 



T T' 







T t T' t 



The 



respectively, which satisfy a^, a^. ' = dxT'^ik — k') and 
field rescaling a~'^ in equation @ can be anticipated as a consequence of broken confor- 
mal invariance. The gauge field quantization (j^) is analogous to the conformal vacuum 
of the scalar field theory in inflation. In the case of a scalar field in conformal space- 
time the procedure is quite straightforward. The kinetic term has a non-canonical form, 
^/~gC]^m{4>) = ci^V^'^ {d^(f))d^(t>) /2. The following simple rescaling dr4> = dr^p/a brings it 
to the canonical form, y/—gC\^in{f) V^'^{9fif){duf)/2. Since the canonical momentum is 
now n^p = dr^p, one can quantize ip by making use of the standard canonical quantization, 
[(f{T,x),TT^{T,x')] = i6{x — x' ). The Drummond-Hathrell action on the other hand 
contains quite a nontrivial modification of the kinetic gauge term. How to determine the 
power ( in equation (||) is discussed in some detail below. 

From the Friedmann equation (|100|) and p oc a^'^^^^'^\ one can infer drd/a = 2/(l-|-3w)r, 
a oc r^/^^+^*"-' (in de Sitter inflation dra/a = — 1/r and a = —1/Ht), so that equation (|) 
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can be recast as 



'^l-{l-6^)V'-'^-^^^AJ = Q. (7) 



r 



where [] 

5^ = 2/?e(10 + 3w7)^ 

e., = - (^^3^1 = -2/3e (26 + 33^/;) ^ (8) 

are slowly varying functions of conformal time. Since oc p oc ']-^^(^+'^)/(^+^'^) ^ this is 
indeed so in slow-roll inflation, for which w ^ —1. When the time dependence of the Hubble 
parameter is neglected, equation (|^ can be reduced to the following Bessel equation 

(9j + ,-^_ML+M)^, = 0, (9) 

where the mode functions are defined in equation (H), and we set 5^, 0. This is le- 
gitimate for superhorizon modes, since the effect induced by 6^ is suppressed by k"^ when 
compared with the ^^„-contribution. In Appendix A we also consider the scalar electrody- 
jl5| and show that for the minimal coupling ^ = and when ,^ < 1/6 the gauge fields 



namics 



get damped in inflation. In the conformally coupled case = 1/6) there is essentially no 
effect, while in the case when ^ > 1/6 the gauge fields grow in inflation (c/. Eq. (|106|) ). 

In deriving equation (|^), the requirement that the coefficient of the damping term drAf 
vanishes leads to the following conformal mode rescaling a~'^ oc (— r)^"" in equation (^, and 
hence 

. (l + 3m) „ 

a-^ = a 2 (10) 



With this conformal rescaling the quantization rules for the photon field in are canonical, 
that is [aT, aV t] = S^^'6{k - k') and the Wronskian W[A'j^\Af] = i. 

The solution to equation (|^) can be conveniently expressed in terms of Hankel functions 
as follows 

= \^r^H^J\-kr), J = 1,2, u = e^ + ^ (11) 



^In fact both 1 — 5w and 9w in Eq. (||) should be muhiphed by (1 — Lw)~^- This however generates 
terms that contribute to the effective action at higher orders in the expansion, so that they can be 

consistently omitted. 
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which, at early times in inflation when r — oo, reduce to the conformal vacuum solutions 

,^oo 1 gTi(fe.+W2W4) + o{{-kT)-^) , J = 1, 2, (12) 

\/2k 



with the standard Wronskian wave function normalization 



W 



A^\Wp]^Wi%A^^-i^.A^^)A^^=^. (13) 



k k \ k / k 



When writing Eq. (pH]) we assumed that Ow > —1/2, so that z/ > 0, which holds in situations 
of physical interest. 

We are primarily interested in computation of magnetic fields on cosmological scales 
today, which correspond to superhorizon scales at late time in inflation, for which fcphys 
For these scales A;|r| ^ 1, and one can use the small argument expansion for the Hankel 
functions -ff^^'^^ in (|n|) to get 



^(1) _ ^(2) * _ 



k k 



.Tiv) f 2Y l-icotTTi// A;rVl ,2.\ 



vr \ kr J r(z/ + 1) 

where we assumed k\T\ <^ 1, u > and z/ 7^ integer. The expansion of Hankel functions for 
u = integer can be found for example in Ref. p3| . 

To model the gauge field evolution in radiation and matter eras, two important effects 
must be taken account of : (a) the energy density p begins scaling away quickly with the 
expansion of the Universe, p oc 1/a^ oc 1/t^, and, (b) as a consequence of the inflaton 
decay into charged particles, the plasma conductivity grows large. Usually the transition to 
radiation era is fast, and typically takes less than one expansion time, while the inflaton decay 
to radiation is model dependent and may take many expansion times, so that one should 
consider these two effects separately. It is important to make this distinction because, as we 
shall see, a consequence of the inflation-radiation matching may be the photon field growth, 
while a large conductivity induces freeze-out of the magnetic field and decay of the electric 
field. 

2.3 Photon field in radiation era 

We shall now compute the gauge field in radiation era in the sudden transition approxima- 
tion. That means that we neglect the effect of conformal invariance breakdown, which may 
play some role at early stages of radiation era. In addition we assume that the effect of 
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conductivity and preheating can be neglected. (For a heuristic account of the conductivity 
at preheating see Ref. |^.) The matching conditions then become simply 

^f{n) = A = «g4+^(r.) + /3!,^4-H^r) 
drAf{n) + ''-^^Af{T,) ^ -^o = «A4^H^r) + /3!,^5.4~^(r,), (15) 

where 4'''^ ~ (2/c)~^/^e^*'^'^ are the asymptotic (conformal) vacuum solutions in radiation 
era. 

Determining the conformal times at the end of inflation (r = Tj) and the beginning of 
radiation era (r = r^) requires some care. We start with the physical requirement, that the 
photon field is in the (Minkowski) vacuum for fcphys = k/a> M. We are free to set a = 1 at 
the conformal time when the vacuum changes from Minkowski to conformal, that is when 
the terms in Eqs. (|I|) and (^ that break conformal invariance become operative. The modes 
k = M then exit the horizon (fc|r| = 1) at a = M/H. This then implies the following 
functional form of the scale factor in inflation 

a(r) = |(-Mr)^/(^+3-) , r<n = -^. (16) 

The smooth matching a(rj) = a{Tr) and dra{Ti) = dra{Tr) then implies for the scale factor 
in radiation 

, , M T l + 3w 1 

4r) = --, r = - — (17) 

The nontrivial term (z/ — l/2)4^''('^«)/Tj in the second matching condition ([T5|) comes 
from the scale prefactor in equation (|^). In this case the coefficients and correspond to 
the Bogoliubov coefficients because we are matching onto the conformal vacuum of radiation 
era, such that the k-mode amplitude amplification can be characterized by the number of 
produced particles Nj: = The result for the matching coefficients is 

Ao-iSo/k Ao + iSo/k 



k 

and hence 

So 

k 

Now taking a small-argument expansion of the Hankel functions (|I^) and assuming u > 
we obtain 



A^{t) = Aq cos k(T — Tr) ^ siu k^T — T^) . (19) 
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Ao = -^r(z/)(7rM)-^(^y , u>0 

/ , 1 /2 (1 — cot nu) ( k 
f„ = i^Mf' ^ ^^^^ > (^—j , 1>.>0 

£o ^ ->r(. + !)(-) (jjj) . .>i. (20) 

where we made use of r(l — 2;)sin7r2; = TT/r{z), such that these equations are vahd for 
integer i/'s as well. To get the photon field amplitude @ in radiation era, multiplication by 
a factor = (H/M) 2~i'^^2) = (^H / M)^^^^^ is required. When expressed in terms of 
Ow = ^ — 1/2, equations ( ^OD become 



j 1 /2M\ t+^^ 1 

A = --r(^^. + l/2)(7rM)-5(^— j , ^»>-2 



1 



1/2 (l^^tanTT^ / ky^ l^fl ^ 1 

M\i/2 / ^ \ l-e™ 1 



£:o = -.r(^. + 3/2)^-j ^— j , ^.>- (21) 

From these and equation (|1^) we easily conclude that the contribution of dominates, and 



hence in this mechanism the electric field contribution £q can be neglected. For this to occur 
the term {v — l/2)A^^\Ti)/Ti in equation (|T3p plays an essential role, since it cancels out the 
leading contribution to £o, making thus the Eq term in equation ([T9|) sub leading. This means 
that all of the amplitude growth by the Drummond-Hathrell action is effected in inflation, 
and the magnetic field strength can be approximated by 

^g(r) = AcosA;(r-r,) (22) 

and is independent of conductivity in radiation era. This is to be contrasted with the 
backreaction mechanism of Ref. |21|], which corresponds to 6'^ + 1/2 ^ M\/H'^, where Ma 



is the photon mass in inflation, such that So/k ^ Aq for superhorizon modes. 
2.4 Magnetic field spectrum 

We shall now compute the magnetic fleld spectrum from the Drummond-Hathrell action in 
inflation and discuss it's cosmological relevance. Since the fleld originates from the amplifled 
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vacuum fluctuations in inflation, it is classical and stochastic on cosmological scales, so that 
the field can be completely specified by the spectrum. By making use of the volume averaging 
procedure described in Appendix B, we find that the gauge field spectrum (^) - ( p2D 
and equations (|T^) and ( |114D imply the following spectrum at the end of inflation 



2 9 X 2^^--^ r (2 - e^) r 



'»+^ 7r3/2 r(^^ + i/2)r(^^^ + 2) ■ ^ ' 

To obtain the magnetic field today, we assume the field to be frozen in the plasma such 
that B(Tq) = B{Th){Th /TqY , and the comoving scale 4 scales as £c = (Th/Tq)^,, where 
To = 2.75 K = 2.37 x 10^^'^ GeV and Th ~ 10^^ GeV is the temperature corresponding to 
the end of inflation when H ^ 10^^ GeV. With this the spectrum today can be written as 

5A.» = -T7^r(^. + 1/2)6,„^J-) -— -. (24) 



When the effects of turbulence are included, a modest amplification of the amplitude on large 



scales may result, provided the spectrum slope 6^ — 2 < —1/2 (c/. figure |^ in Ref. pT|| ). 
To evaluate the magnetic field strength on cosmological scales, the following conversions are 
useful : (lGeV)2 = 1.44 x lO^^ Gauss and lOkpc = 1.56 x 10^^ GeV~^ In figure | we show 
the spectra obtained for various 9yj of physical interest. In particular, the vacuum spectrum, 
which corresponds to the conformally invariant case, is recovered when 6'^ = 0, and can be 
inferred from equation ( |116D : 

B}^ = J- £-2 ~ 3 X 10-55 (^^^^^ Gauss (vacuum) . (25) 

The thermal spectrum corresponds to = 1/2 in equation ( ^4]) and reads 



4(1+-) ■ • - 



Since the lower bound on the seed field for the galactic dynamo is about i?sccd > 10 '^^ Gauss 



at a comoving scale of £c ~ 10 kpc [^, |2^, both of these spectra are too weak to be of 
cosmological interest. When 6^ = ^ however, equation (0) reduces to 



HJ 
11 



which is sufficiently strong to seed the galactic dynamo. In other words, to obtain the field 
of cosmo logical interest, the spectral index must be greater than about —1. An exception 
are models in which the photon becomes massive dynamically in infiation [^, Finally, 
when = 2 the spectrum becomes scale invariant with the field strength (c/. Eq. ( |1 1 7| ) ) : 



E;;-^""- ^ 3 X 10"^^ l^—j Gauss . (28) 

Here we have not displayed the logarithmic enhancement ~ (— ln2A;o^)^ in Eq. (|117|) , where 
ko is the smallest momentum amplified in infiation. Note that for M H the field (^) is 
somewhat weaker than B ~ 10^^ Gauss suggested in 0. Furthermore, it conforms with the 
magnetic field bound from the CMBR considerations in [^]. In figure ^ we illustrate all of 



the spectra (|25|) 
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Figure 3: Magnetic-field spectra (p^) from the Drummond-Hathrell action as a function of the 
comoving scale £c today. The vacuum spectrum (|25| ) {6w = 0) is shown in dot- dot- dashed green, the 
thermal spectrum (|2^ ) {6^ = 1/2) in dot-dashed gold, the 9^ = 1 spectrum ( |2^ ) in dashed red, and 
the scale invariant spectrum ( pSj ) {9^ = 2) in solid blue, where we assumed that the fermion mass 
M H. We also show the relevant dynamo bounds -Bseed ~ 2 x 10^^'' Gauss for a universe with 
critical matter density, and -Bseed ~ 2 x 10~^^ Gauss for a flat, low-density universe with dark energy 
component. Note that when 9^ ^ 1 the dynamo bound -Bseed ~ 2 x 10~^^ Gauss at ic ~ 10 kpc in 
a flat, low-density universe is amply satisfied. 
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An important question is of course whether one can obtain 6^, ^ 1 from the Drummond- 
Hathrell action with a reahstic choice of parameters. To investigate this, we write equa- 
tion (IP as 

(A±J^^imj^?!^^^e.. (29) 

[1 + w) ae rl^ 

To good approximation, 6^ in equation (|^) maximizes when w ~ = —59/66 (for which 
— (26 + 33w)(l + w) = maximum), such that 

, 16657r , , 

Np{B^ = 1 ~ — , 30 

and hence ~ 10^ charged fermions with a mass M ^ Ji resuh in 6'^ ~ 1, too high to 
be expected in grand- unified theories. In this case lyj ~ 30, indicating break-down of the 
expansion. Further, in this mechanism the value = 1 is not in any sense favoured, 
or naturaL In conclusion, we have shown that is it quite unlikely that magnetic fields 
can be amplified by the Drummond-Hathrell action in inflation to a strength of interest for 
cosmology. Moreover, our result should be interpreted with caution, since large amplification 
is obtained in the region of parameter space when the expansion becomes unreliable. In 

order to address this question one would have to study systematically higher order corrections 
in the Schwinger-DeWitt expansion |1H . 



3 Photon coupling to gravity with light fermions 
3.1 Trace anomaly 

Dolgov 1 10] has argued that conformal invariance of gauge fields gets broken by the one-loop 
diagram shown in figure |1|, even in presence of massless fermions. Taking account of the 
local contribution only, in a conformal space-time with the metric g^^y = a^rifj,^, the following 
effective action 

5anom = J d'^x{-g)^^'^ 



(31) 



where 



i^eiNp) = — (32) 

ovr 



reproduces the equation of motion in Ref. |[T0[- This action manifestly breaks conformal 



invariance, since it is not invariant under the transformation a Xa. It is however invariant 
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under the following combination of conformal transformation and dilatation, a Aa, and 
x^" A'x^, with A' = X~^-/^. 

When non-Abelian gauge fields with an SU(N) symmetry are considered, the anomaly can 
be studied by {Np, N) = 2a(.{Np — llA^/2)/37r. In the limit when the fermion mass 

<^ |7^|, where TZ is the curvature scalar, one can partially resum |Q terms involving the 
curvature scalar in the Schwinger-DeWitt series |^5|, |15|. For example, the effective action 
of a non-Abelian gauge theory contains the following logarithmic term 



Sin = J d^x{-g)2 



11 a,C, ( 7^ 



-g^^'^gP'^F^^Fl 1 ^ In , 



(33) 



where /i is a mass scale, = e^/An, F^^ is the gauge field strength, and C is the defined 
in terms of the structure constants by trlT^'T^] = -f^^^jbcd ^ _CS''\ such that C = for 
SU{N). When the curvature in (^) is large, the effective gauge coupling becomes small, 
and one encounters a curvature induced asymptotic freedom. Now, since in inflation TZ = 
-(1 - 3w)p/M^ {cf. Eq. (|9D) and p oc a-^^^+'"\ we have TZ-^dn/dr = -3(1 + w)da/adT. 
This then implies that one can include the effect of the logarithmic curvature running from 
the non-Abelian loop corrections in inflation by the replacement Kg — ^ K,e + i^+w){llaeC/A7!-) 
in Eq. (|3^ {cf. Eq. (|3^)). Similarly, for scalar electrodynamics we have Kg ^ /tg + (1 + 
w){aeNs/87i), where A^^ is the number of light scalar particles (m^ <^ l"^!)- With these 
remarks, the analysis presented in this section applies also to the resummed scalar curvature 
corrections. Note that in de Sitter inflation w = —1, and hence the logarithmic corrections 
in (|53| ) do not lead to any amplification of gauge fields in inflation. 

We will now show that, just as in the case of the Drummond-Hathrell action, in inflation 
the gauge fields are amplified in presence of massless fermions, and damped in presence of 
massless bosons. Making use of the variational principle, SSi^riom/^^u = 0, we arrive at the 
following equation of motion for the transverse gauge field 



d''^-V' + Ke—dr]A^ = 0. (34) 



a 



We now quantize the field as in equation (|^) and get the following mode equation 

«^-*-=-f(-?)(^y-f?)A-0 (35) 
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with the conformal rescahng 

Assuming a power law inflation with the scale factor a oc r2/(i+3«') (see Eq. (|16D), this 
becomes 

{dl + k'-^^^^^]A, = ^, < = -^, (37) 



such that Kg can be identifled with 9^ of equation (j^), which is also consistent with equa- 
tion (plQl). This than means that the mode functions in inflation are solved by equation ([TT|) : 

Af = \^^Hil\-kT), J = 1,2, v, = \ + k',, (38) 

which satisfy the Wronskian (plsl) . This is consistent with the canonical quantization: 

[Ai{T, x) ,Iij{T,y)] = i6ij6{x — y), where Uj = a'^^drAj denotes the canonical momentum. 

3.2 Photon in radiation era 

An important difference between the Dolgov anomaly and the Drummond-Hathrell action is 
that the anomaly also survives in radiation and matter eras, since it is does not dependent on 
the energy-density scaling. As a consequence the matching conditions on the mode functions 
now differ from Eqs. (|15]) : 



A'V{n) ^ A, = a^,A'p{r,)+(3*_^,A\'yn 

drAf{T:^ = -So = aj^drA^^\T,)+Py^d^At\Tr), (39) 



where the radiation era modes are the following spherical Bessel functions 

k ^ ' k ^ ' 2 Ur \ / v2A; 2 



The following parametrization of the scale factor is convenient (c/. Eqs. (|T6|) - (jl^)) 

We now make use of the Wronskian W[^g^'', ^i""*] = i, and arrive at 

c^k = -^[Af{n)^^A7\r^)-{^U^{n))A^-\r^)_ 

= ^[Af{T.)d^A^^\rr)-{d^A'i\n))A'^^\T^)] (42) 
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Vi-Vr 



which, with the help of equation (0), reduce to 

+ o(r>-^'-)- (43) 

For Kg > this is the dominant term since in this case i/j — i/^ = — (1 — 3w)Ke/2(l + 3w) > 0. 
Note that amphfication maximizes in de Sitter inflation, when = —(3_j: oc kr^'^. After 
some algebra we obtain for the gauge field in radiation 



-AfM = 4'"' COS (^tT - ^(^, - 1/2: 



TT 



where we used equations ( ^0|) and (|43|) . For de Sitter inflation this reduces to 



^r(^) ^ ^ ^ ^ cos(fcr + 7:kJA) . (45) 



^/2k vr V A; / 

When conductivity grows quickly in radiation era to a large value, the smooth matching ( ^91) 
is not any more appropriate. In this case the field amplitude freezes-out to the value at the 



end of inflation (|38|) implying that, at the beginning of radiation era. 



Note that this amplitude is smaller than (^41), since {vi — Vr) — ('^i — 1/2) = Ke/2 > 0. 
3.3 Magnetic field spectrum from the Dolgov anomaly 



We first consider the simpler case (^Bf), for which the magnetic field gets frozen in the plasma 
at the beginning of radiation era. The volume-averaged spectrum can be computed as in 
section |2.4|, resulting in the following spectrum today 



2 



^ 9 X 2-^-^ r (2 - kJ2) V {kJ2) 
^ 7r3/2 r(K + l)/2)r(2 + /te/2)' ^ ' 

where for simplicity we assumed de Sitter inflation, for which the resulting field is strongest. 
The scaling (p6D introduces ambiguity to the field amplitude normalization. Our normal- 
ization corresponds to setting the scale factor a{Ti) = 1 at the end of inflation, just as in 
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equation (0), such that at r = Xj the horizon scale modes have the amphtude of the order 
the standard Minkowski vacuum amphtude. We beheve that this justifies our normahzation 
prescription. The analysis of the spectrum { ^7\ ) is identical to that of the Drummond- 
Hathrell spectrum (^) discussed in equations (^) - with the replacement 6^ —>■ i^e/'^ 
and M = H. 

The second case (^) - (145|) applies when the conductivity is low in 'radiation era'. This 
means that the infiaton oscillates for a long time before charged particles are produced that 
freeze-out the magnetic field. Denoting the freeze-out temperature by Tf^eeze T^, we 
then find that the gauge field dynamics in de Sitter infiation ( ^51) induces the following the 
spectrum today 

^ 9 X 2^--« r (2 - /.e) r (/.e) , . 



There is subtle difference between this spectrum and (p4D. The additional suppression here is 
due to the anomalous scaling a^^l^ = {Tireeze/TnY^^'^ of the gauge field. Here Tfreeze denotes 
the temperature at which the modes freeze out due to a large conductivity. This expression 
is correct for the modes that enter horizon before the freeze-out. We do not discuss here the 
modes that grow for a while and freeze out before reentering horizon in radiation era. The 
dynamics of these modes is quite complex and their spectra lie in between those given in 
equations (g^ and (|||). 

The physics of the infiaton decay may in fact be nonperturbative, and there is a large 
literature on the infiaton decay and preheating after infiation. Nevertheless, the gauge field 



dynamics at preheating is not yet fully understood [Q. A heuristic account of preheating 
with gauge fields is given in Ref. 



The magnetic field strength corresponding to the spectra (|47|) and is very similar 
to those plotted in figure ^ and hence we do not plot them here. An interesting question 
is of course under what conditions p7|) and ( ^HD can result in spectra with a slope >, — 1, 
required to seed the galactic dynamo. To get the spectrum Bi oc equations ( ^71) and ( ^H]) 
imply 1 < Kg < 2, or equivalently 

^<Ar^<^. (49) 
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For ae ~ 1/40 this gives the following bound 

2 X 10^ < AT^ < 4 X 10^ 



(50) 



for the number of charged fermions. In presence of bosonic fields amplification is reduced. For 
example, when the electromagnetism is embedded in an SU(N), then equation (^9]) should 
be read as a limit on Np — UN/ 2. When applied to the standard model, the photon field is 
amplified on scales smaller than the electroweak scale My/ ~ 80 GeV. On higher scales, the 
hypercharge field is amplified, and the relevant number of light fermions is Np = 12 (above 
the top mass). Hence this mechanism results in weak magnetic field amplification if the 
standard model fermions are the only light fermions. To get strong amplification, one would 
need a few hundred of additional fermions with masses at some intermediate scale below 
inflation. It is not clear whether this can be made consistent with the grand unification. 

The validity of the Dolgov anomaly equation ( |5^ ) was questioned in [|T^. Here we have 
solved Eq. (^) exactly to order e^. To properly address the question of validity of this 
approach one would have to study the higher order corrections in the coupling constant 
expansion to the trace anomaly. 



4 Photon coupling to scalar and pseudoscalar fields 
4.1 Photon coupling to scalar field 

Consider now the following action with gauge invariant coupling of the gauge field to a scalar 
field 

S^ = j dS^g ^-^-f{^)g>^Pg--F,,F,^ + ^g'^''{d,<P){d.<P) + ^(0)) , (51) 

where F^j, is the gauge field strength, is a scalar field and V{(j)) is the corresponding 
potential. In general (p may couple to other fields as well. The 'coupling' function / = /(</>) 
we take to be of the form: 

where /„ are some coupling coefficients and Mp = {8ttG)~^^'^ ~ 2.4 x lO^^GeV is the reduced 
Planck mass. For the standard gauge kinetic term, /o = 1, and in addition the potential 
V{(f)) must be chosen such that today = and V{(j) = 0) = 0. A particular realization 
of this model has been considered by Ratra [|lT|, where /(0) oc e"^'^^^^'', and was taken 
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to be the inflaton in an extended inflation with the potential oc e ^'i'l^^p _ In this 

case fn = (— A)"/?T,!. Ratra's model has subsequently been reconsidered in the context of 



string-inspired cosmology []T2|, where the role of is taken by the dilaton field. Here we 
study the conditions under which /(0) and V{(f)) can lead to significant amplification of the 
magnetic field in inflation. As a particular case we discuss Ratra's model. 

Working in Coulomb gauge, Aq = and V ■ A = {cf. Eqs. (H-H)), the action implies 
the following equation of motion in de Sitter inflation for the transverse gauge field 

{d'r - V' + ^ (5^1n/(0)) d)j Aj = 0, (53) 

where we assumed that in inflation the field (p = (/)(r) is a function of time only. It is 
now clear that the coupling function / is responsible for breakdown of conformal invariance, 
and hence its precise form determines whether the field gets amplified in inflation. In what 
follows we consider two types of functions: (A) quadratic coupling function 

and (B) exponential coupling function 



/ = l + /2^ + - (Type A), (54) 



/ = e-^*/^- (TypeB). (55) 

Type A coupling functions are expected to be induced by the Planck physics. The coupling /i 
can be made to vanish in (Q) by imposing the symmetry —0. Similarly, by promoting 
to a complex field and imposing a Zn symmetry, one gets that /„ = 0, Vn < A^. Type B 
coupling functions may occur from string theory corrections to the gauge field dynamics. 

4.1.1 Chaotic inflation 

Making use of equation ( p.23|) in Appendix C.l where we estimate drcf) during slow-roll in 
chaotic inflation, where V = m^0^/2, we can recast equation (^3]) as 



T 

For Type A coupling functions this implies 



_ + (a,ln/(0)) Or] Aj = 0. (56) 



- - Af = 0, 9, = -2/2, (57) 



r 
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where we neglected the higher order corrections. Similarly, for Type B coupling functions we 
obtain the following equation: 

[d^r - - ^5.) A- = 0, e, = A^, (58) 

where = 0o + (2M|,/0) \n{—HQT). Since 62 and 9\ are generally slowly varying functions 
of conformal time, equations (|57| - ^8D can be well approximated by the Bessel equation (0) 



analysed in detail in section |2.1| . The corresponding magnetic field spectrum is calculated 
in (p3|-p^, according to which the spectrum reads 



E,, cxC+e^ {6 = 02,9^), (59) 

where is the comoving scale. To get the amplitude one should set M = if in (^3|-p^), 
which corresponds to the normalization to the Minkowski vacuum at the horizon crossing. 
We assume that this normalization is effected by the local interactions in the de Sitter vacuum 
operative on subhorizon scales. With this the magnetic field spectrum can be simply read 
off from figure 3. To get an acceptable spectrum, we then have 2 > 62, 0\ > 1, or 

\< -f2 < 1, A - 30, (60) 

where we took Mp/0 ~ 1/(2 x Q^/'^N'^/^) ~ 1/20, and the number of e-folds ~ 50. We 
have thus found that one can obtain magnetic field spectra with observable consequences 
with /2 of order unity, which is a natural value. A generalization of Type A model contains 
higher order terms in the 0/Mp-expansion. Since in chaotic inflation typically cf) ^ Mp, 
the resulting magnetic field spectra of interest can be obtained with coupling coefficients 
|/2n| <^ 1 (?T. = 2,3, ..), where /o = 1, /i = /2 = /s = 0. On the other hand, for Type B 
model the coupling A is required to be unnaturally strong. Moreover, one gets a significant 
amplification in inflation only when the effective gauge coupling e^(0) = Cq e"'^'^/^^^ starts 
very weak in inflation. 

In the spectrum computation in section ^.4| we have neglected the conductivity in radia- 



tion era. A simple way of accounting for that effect is to assume that magnetic field freezes 
in when the conductivity becomes large. The spectrum calculation in this case can be per- 
formed analogously as in section {cf. Eq. (|3))- To get spectra of physical interest in this 



case, one requires by about a factor of two larger values of 62 and 9\ from those indicated 
in (H). 
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4.1.2 Extended inflation 

For the potential of extended inflation V = Vqc"^''^/*^^ equation (|53D can be written as 

dl-V'- y^^^l id,\nm) d}j Af = 0. (61) 

where, to estimate d(f)/dT, we made use of equation ( |131D in Appendix C.2. For Type A 
model dU) this implies 

where is a slowly varying function of conformal time given in ( |132| ). Since ~ 0o ~ 
— {Mp/X') ln(12/A''^), this immediately implies 

/ 12 

l<e',<2 ^ /^^-(^in^j (A'«l), (63) 



which is to be compared with the conditions in chaotic inflation (^) and (|60|). For Type B 
model (|55|) we get the following equation: 

To get a magnetic fleld of observable strength we than have 

1<^1<2 =^ ^<-A<^ (0<A'<1). (65) 

In this case the effective gauge coupling e^(0) = e^e~^'^/^'^ begins very weak at early stages 
in inflation (0 ^ —Mp) and acquires today's value at the end of inflation (0^0). Precisely 
this feature of the pre-Big Bang model was used in [1^ to get magnetic fleld ampliflcation, 
where the gauge coupling constant is determined dynamically by the dilaton expectation 
value. 

An interesting question is whether there is a natural spectrum predicted by inflation. 
In other words, is there a value of the spectrum slope 2 — 9 that is in any way singled 
out. Assume that at early stages of inflation 6 > 2, which implies a (tachyonic) instability 
resulting in a fast amplitude growth and spontaneous magnetization on superhorizon scales. 
This condensate may affect evolution of the flelds in inflation and possibly single out the 
scale invariant spectrum. This question deserves further investigation. 



21 



4.2 Photon coupling to pseudoscalar field 

We now consider the following coupling of a pseudoscalar field the Chern-Pontryagin density 

(66) 



where 



is a coupling function of a pseudoscalar field ip, F^^ = efj,i^ripFrip/2 is the dual field strength 
and e^i,rfp is the totally antisymmetric tensor defined by eoi23 = 1- The case h{ip) = h2ip^/Mp 
is mentioned in ^ as potentially interesting, but no attempt was made to solve it. There 
is an observational constraint on h coming from the magnetic dipole moment measurements 
of the neutron and electron, according to which |/;-o| ~ 10"^ [ip = 0); for simplicity here 
we assume Hq = and ip = 0. With ip = today there are no experimental constraints 
on hn {n > 1). When applied to the quantum chromodynamics, the question of smallness 
of h is termed the strong CP-problem. Peccei and Quinn have suggested that h might be 
dynamically driven to zero by symmetry breaking of ip at some high energy scale. The 



special case of this model was considered in |T^, where the authors assumed that the axionic 



potential is of the form V{tp) = Vo[l — cos^ip/ipo)], and found no significant amplification. 

The action ( |5BD implies the following equation of motion for the gauge field in inflation 
in Coulomb gauge (c/. Eq. (^)): 

{dr - V' + ^ In h{^))Vx^ Af = 0, (68) 

where we used rj^^'d^Fpi, = and assumed that in inflation ip = iP{t). Conformal invariance is 
now broken hy h = h{tlj), which then may lead to magnetic field amplification in inflation. It 



now follows from the analysis in section [4.1| , that in inflation quite generically equation ( |68D 
can be rewritten as 

(^d'^^V'-^Vx^A^ = 0, (69) 

where O^p is a slowly varying function of time in inflation, and 6'^ ^ in radiation era. The 
specific forms of 6'^ for Type A and B models (5^-^ in chaotic and extended inflation can 



be easily reconstructed from Eqs. (|57| - p8D , (|6^) and (|6^) , provided f h and (j) —> ijj 
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Equation (^) can be written in terms of the mode functions (|^) as follows 

(a? + *^)44-T^*x44 = o- (70) 

Since the circular polarization vectors satisfy ik x e"^ = ±ei (T = ±) we can write (|7Dp for 
the circularly polarized mode amplitudes as 

[^9^ + P T '^k^ 4 = 0. (71) 

Note that when < {6^ > 0) the mode A'^ (-^^^ S^^^ amplified in inflation (where 
we took into account that r < in inflation), such that in principle one can distinguish 
the magnetic fields produced by a pseudoscalar coupling to the Chern-Pontryagin density 
from other mechanisms discussed in this paper, where the field amplification is polarization 
independent. This is of course true provided amplification is strong enough to be observable. 
When written in terms of the variable u = 2ikT, equation (fflD reduces to the following 



Whittaker equation 33 



(aj_l±!^)^. = 0. (72) 

The solutions can be expressed in terms of the Whittaker functions Wi ;^{u) as 



with the Wronskian normalization W[^^ (r), ^5 (— r)] = i, where ^5(— r) denotes the (sec- 
ond) linearly independent solution. The asymptotic form of ([73|) is 



^ \/2k 



such that it approaches the Minkowski vacuum at r ^ — oo. The only effect of the ■j/'-field 
at a distant past is a time-dependent phase shift in the mode functions. For superhorizon 
modes close to the end of inflation the following asymptotic form is useful 



Now matching this smoothly onto the radiation vacuum solutions A^'''^ = (2A;)-i/2eTifcr 
results in the vacuum spectrum of a modified amplitude. The circularly polarized modes A"^ 
get amplified in inflation by a factor 

A 



k 



itrad 



elie.lTf^v-, (76) 
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where we made use of the Stirhng formula, |r(l ±i6'^)| ~ (27r|6'^|)^/^e~''l^"^l/^ {\6^\ > 1). 
This means that, when 9^ < { 9^ > 0), the mode i^-j:) gets amplified by a factor 
~ g^ie^l^ while the amplitude of A'j^ i-^^) remains unchanged. This is an extreme case of 
birefringence, where the photon of one polarization corresponds to a tachyonic particle, and 
the other to a massive vector particle. 

We have thus shown that a pseudoscalar field coupling to the Chern-Pontryagin den- 
sity (|66|) cannot change the vacuum spectrum Bi oc of gauge fields in inflation. In this 
case conformal invariance is broken in inflation such that the vacuum amplitude changes 
and an unobservable phase shift to the gauge mode functions is induced. The resulting am- 
plification of the vacuum spectrum is typically too weak to be observable today. Moreover, 
thermal processes in radiation era are expected to create magnetic field with the spectrum 
Bf^ (X [^, which dominates over the vacuum spectrum on scales of interest in cos- 

mology. 

Even though this model is not relevant for cosmological magnetic fields, it may be of 
interest for baryogenesis. An explicit CP-violation in the pseudoscalar sector may result in 
a nonperturbative amplification of the CP- violating phase of a complex pseudoscalar field at 
preheating. When the field decays into the standard model fermions, a net baryon number 
may be produced. This is reminiscent of the Affleck-Dine mechanism [^, with an important 



difference however: the ip-field is a pseudoscalar field which is not charged under the baryon 
number, so that baryon violating processes are needed to get a nonzero baryon production. 
In a related work the dynamics of a pseudoscalar field in radiation era and its relevance for 
baryogenesis has been considered in 121, |3( 



5 Photon coupling to scalar cosmological perturbations 



An interesting proposal has been recently put forward by Maroto pO[, who considered the 
effect of the photon field coupling to the scalar metric perturbations. We now reconsider this 
mechanism (without making a perturbative expansion of the mode amplitudes). In presence 
of the scalar cosmological perturbations the metric gets modified as follows 

g'^'^ = ((1 - 2<l>)rfr2 - (1 + 2^)dx^) , (77) 
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where $ = $(77, x) denotes the (gauge invariant) scalar gravitational potential. The potential 
$ is to good approximation scale invariant. The CMBR measurements and the large scale 
structure formation imply that $ ~ 10~^. 

The gauge field equation of motion is simply 

d, {V^g^'g^^F,^) = 0, (78) 

which, when written in components and working to linear order in $ (|$| <^ 1), reads 

d,[{l-2<l>){diAo-drA,)] = Q (79) 

and 

dr [(1 - 2^){diAo - drA,)] + dj [(1 + 2^){d,A, - diAj)] = 0. (80) 
Equation ( |79D can be also written as 

VMo - 2V$ ■ VAo = drV-A- 2(V<I>) ■ drA (81) 

where we used Vln(l — 2$) ^ — 2V$. Setting Aq = leads to the following generalized 
Coulomb gauge 

drV ■ A = 2{V^) ■ drA. (82) 

It turns out that this modification becomes unimportant on superhorizon scales, where 
fc|r| <^ 1. We can now write equation ( pO|) in the gauge (|8^) as 

[dl - V') Ai - 2{dr^)drA, + d^V-A- 2{dj^){djA, - d^Aj) = 0, (83) 

where we neglected $ in favour of one in 1 — 2$ and 1 + 2$, which is a controlled truncation 
when working to leading order in $. The last two terms in (^) can be consistently neglected, 
since formally they are a small correction to the V^Aj-term. This follows from equation ( p2D 
and $ ^ 1. Moreover, it is easy to see that the term 2{dr^)drAi (which potentially leads 
to gauge field amplification) dominates on superhorizon scales {k\T\ ^ 1) over the terms we 
have just argued can be neglected. With this equation (|83| ) simplifies to 

- V^) A, - 2{dr^)drA, = 0. (84) 

We shall now solve this equation in the mean field approximation 

$ ^ $^ = ($)^, (85) 
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where ($) l denotes a spatial average over some scale L. Since $ is almost scale invariant 
in inflation, we expect $i to be a slow (logarithmic) function of time and L. Taking for 
example ~ <l>oln(— r/L) ($o ~ 10"^ ^ 1), we find 9,-$ ^q/t. With this we can write 
equation ( |8^ ) in the mean field approximation as the following Bessel equation 

dl-W^- ^9.) A, = 0. (86) 



In section ^]T] we study in detail this equation (c/. Eq. (|^)) and show that the solutions 



for the mode functions can be written in terms of the Hankel functions (ITTl) with the index 



I + $07 and that the resulting mode amplification oc {—krY is small when 



\v — 1/2| = |$o| ^ 1- Since in radiation era 9^$ = 0, the mode functions are simply 
A^^^ = (2A;)^^/^e^*^'^, so that the calculation in section |2]^ leading to the spectrum ( p3D 



applies, and hence we get an almost vacuum spectrum today Bg oc ($q ~ 10~^). 

We shall now argue that the mean field approximation represents an upper bound on the 
amplification on superhorizon scales. This approximation captures the effect of a spatially 
averaged gravitational potential, but does not account for inelastic scatterings. Since inelastic 
scatterings are in general dissipative, their inclusion can lead only to suppression of gauge 
fields. Moreover, the sign of in (|53| ) can be both positive and negative, suggesting an 
additional suppression when compared to (^^- To show that this is indeed so, we now study 
Eq. ( ^4|) by making an Ansatz for the momentum dependence of the mode functions. Upon 



performing a Fourier transform on (|j) we get 

{dl + A- 2 / ^(a.$,.oa.A^_g, = 0, (87) 

where we suppressed the polarization indices. The mean field approximation consists of 
neglecting inelastic scatterings that lead to nontrivial momentum exchange, i.e. it amounts 

The spectrum of the potential $ is typically close to scale invariant, which means that 

= 0o(r)A^-'fc-^, A; < A, (88) 
where (3 is the spectral index for superhorizon modes (fc < A), and 

<Po = Co-^(- fa'dr) (89) 
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is a slowly varying function of r, and Cq ~ 10^^ is a normalization constant. For a scale 
invariant spectrum (3 = 5/2 — n = 3/2. On large cosmological scales the CMBR measure- 
ments constrain \(3 — 3/2| ^ 0.2; on smaller scales a similar constraint results from large 
scale structure formation. The cutoff A in (RSI) is of the order the horizon H at the end of 
inflation, such that on scales /Cphys = k/a > A H the vacuum spectrum |$^| ~ |0o|/v^ is 
recovered. Similarly for the gauge field we take 



A,: = A(r)A;-°, A; < A, 



(90) 



where the scale invariant spectrum for the magnetic field corresponds to a = 5/2, and 
Aj: oc A;^^/^ for k > aA. Now inserting the Ansdtze (|88| - pOD into (|87D results in 



k 



3-/3 



27r2(2-/3) \A 



— Iaf3iXo, Xi) dr Aq 



where 



XI 



dxx 



l-a 



[l + X 



,2-/3 



— |1 — a; 



2-/3 



(91) 



(92) 



where Xi = aA/k, and we have introduced a lower limit of integration xq = k^/k, where 
ko denotes the lowest moment amplified at the beginning of inflation. The self-consistent 




Figure 4: The itns/iadec? region corresponds to the phase space where Xq,/3(0, xi = A/fc) is dominated 
by the ultraviolet cutoff, and where the only self-consistent solutions of (^ij) lie (marked by the 
thin vertical line a ~ 1/2, |/3 — 3/2| < 0.2). In the light blue shaded triangle the integral Xq,^(0, oo) 
in ( |92[ ) is convergent. 



solution of (|91|) is obtained when when a + P < 3. In this case the ultraviolet cutoff Xi 
^For the special case /? = 2 the integral Xafiixo, xi)/{2 — (3) ^ dxx^^" In [(1 + x)/\l — 
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dominates the integral (p^), and equation ( pT] ) becomes 



Since in inflation d(f)o/dT ~ Cq/t, and a ^ —1/Ht, Eq. (p3D can be simplified to 

(a? + ,^ + 2,.m)A.^0, -,,(,) (94) 

When a + /? 2, 7,- ~ 7o is almost time independent, and equation (0) is self-consistently 
solved by 

I vifc^ '^^'^ when 'Jo > k 



where cuq = yfc^ — 7o and 7± = 7o± y7Q — A;^ > 0. The critical case kc = 7o(fcc) corresponds 
to kc ~ (CQ/47r'^)(_f/'^/A) ^ In all cases in ( P^D the spectrum on large superhorizon 
scales is to good approximation the Minkowski vacuum, and the gauge field actually decays 
slowly in infiation. This suggests that the mean field equation (^) can only overestimate 
amplification, and in addition it indicates that the vacuum a = 1/2 approximates well the 



actual solution. The solution to (|9^ ) is represented by the thin vertical line in figure ^ 
and it is the only self-consistent solution of (0). Indeed, when a > 3 the integral (|9^) is 
dominated by the infrared cutoff. In this case la/sixo, oo)/{2 — /3) ^ (2/(q; — 3))(/c/A;o)°~'^5 
which cannot be made consistent with (^7]). Finally, when a, /3 < 3 and > 3 the integral 
Ta/3(0, oo) converges (in the light blue shaded triangle in figure ^) and can be evaluated |^ 



to yield X„/3(0, cx)) = 5(2 - a, a + /5 - 4) - 5(3 - /3, 2 - a) - B{a + /3 - 4, 3 - where 
5(/i, v) = r{fi)r{u)/r{fi + u) denotes the beta function. Now since |/3-3/2| < 0.2, Eq. 
can yield only a tiny amplification of the vacuum spectrum, such that a ~ 1/2, inconsistent 
with a + P > 3. 

To conclude, in this section we have reconsidered the coupling of gauge fields to the scalar 
metric perturbations and showed that there is no significant amplification of gauge fields 
in infiation. We first analysed the problem in the mean field approximation (which overesti- 
mates the effect), and then we estimated the effect of dissipation from inelastic scatterings 
by employing a simple approximation for the spectra of $ and Ai. The disagreement with 
Ref. |2^ is probably due to the use of the time dependent perturbation theory in [^], which 
often leads to spurious secular effects when applied in large time limit. 
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6 Photon at preheating 



Since at the moment there is no satisfactory account of both the resonant inflaton decay 
and plasma conductivity at preheating, in the present work most of the times we make the 
conservative assumption that no amphfication occurs at preheating epoch. For an heuristic 
account of this problem see for example Ref. and for a numerical work see Since 



explosive particle production is quite typical for many models of preheating, abundant pro- 
duction of charged particles occurs quite generically immediately after inflation. The sudden 
matching approximation |jll| may thus be appropriate, according to which the conductivity 
becomes large quickly at the beginning of radiation epoch such that the electric field gener- 
ated in inflation decays, while the magnetic field gets frozen in. In this work we adopt this 
simple approximation. For a recent reformulation of the matching conditions appropriate 
for treatment of superhorizon charged scalar fluctuations created in inflation in terms 



of more sophisticated fleld theoretical techniques see Ref. |18 



7 Discussion and concluding remarks 

We have considered various effects that break conformal invariance of gauge flelds in inflation. 
We flrst analyse the Drummond-Hathrell action in section ^ which describes a local one-loop 
modification to the gauge field action in presence of massive fermions, and find only a very 
modest amplification of gauge fields. On the other hand in section ^ we consider the Dolgov 
trace anomaly of gauge fields in presence of massless fermions and find a significant gauge 
field amplification in infiation, provided the number of light fermions is of the order a few 
hundred. We have then in section considered modified kinetic gauge terms by coupling 
to a scalar field and found that under quite reasonable conditions on the couplings one 
gets significant amplification of gauge fields in infiation, leading to potentially observable 
magnetic fields on cosmological scales. In section |4.2| we consider coupling of a pseudoscalar 
field to the Chern-Pontryagin density, and find only a mild amplification of the gauge field 
amplitude in infiation, while the vacuum spectrum remains unchanged. We find in this case 
a nice example of birefringence induced by gravity, where one photon polarization exhibits 
tachyonic growth in infiation, while the other, which corresponds to a massive excitation, 
remains of a constant amplitude. Finally, in section | we consider gauge field coupling 
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to scalar metric perturbations and find that the effect is by far too weak to produce 
observable magnetic fields. 



When this work was nearing completion, Ref. |^ appeared which overlaps in part with 



section [4.1| . The results of section [4.1| were presented earlier at the workshop Beyond the 



Standard Model in Bad Honnef (March 2001). 
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Appendix A 

Here we present derivation of the photon field equation of motion starting with the Drummond- 
Hathrell action (0). In a conformally flat space-time the equation of motion can be recast 
as 

V^'d.F,, + -^d, [br^^^'TZF,, + ^a' {TZ^" + Fp, + da^^TZ^^'^ F,J^ = . (96) 

The metric of a conformally fiat space-time can be in general written as g^u = o^{ri)ri^y, 
where r]^y = diag(l, —1, —1, —1) is the Minkowski metric, and the scale factor a = a(r) is 
a function of conformal time r. The nonvanishing components of the Riemann tensor have 
the following simple form 

V ^ — —V ^ — 



7? - _7? 




n,,,' = {v.%k-v^kvm^] (97) 



such that the curvature tensor components are 



TZqq — T^^oio — -3 
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-1->0 _|_ -T->« 

iOj ^ ilj 




and the Ricci curvature scalar reads 



n = -Q 



dla 



-MGp (1 - 3w) 



(9^ 



(99) 



p/p (—1 < w < \). To get the last equality we used Einstein's equation 
.^y lZgfj,y/2 = SnGT^u which, in a conformal space time with the metric g^i, = a?ri^y, read 



where w 
R 



dfa _ dla 



drOX ^ 


8ttG 




3 P 


dra \ ^ 


AtxG 




3 



(Friedmann) 
(p + 3p) , 



(100) 



where we used the ideal fluid stress energy tensor = diag(p, — p, — — p), where p and 
p denote the (total) energy density and pressure, respectively. When p + 3]9<0(— 1< 
w < —1/3) the Universe undergoes an accelerating expansion (inflation); the limiting case 
w = —1 corresponds to de Sitter space-time, while matter and radiation era are obtained 
when w = and w = 1/3. More generally, p oc a~^^^~^'^\ such that upon integrating the 
Friedmann equation one obtains a oc r^/(^+'^"') (in de Sitter inflation a = —l/Hr). 

Now making use of equations (^ - ^) and working in Coulomb gauge in which Ao = 
= V • A, we simplify equation (RBI) to get 



+ 



{dl 

M2 



(66 + c)^ + (c + 2c/) 



dr 



(66 + 3c + 2d)^drA^^ + (3c + 2d) 



( 



2 aT 



0, (101) 



where we used F^i = drA^'^ and djFji = V^A^^, and Aj'^ = (5ij — didj/'V'^)Aj denotes the 
transverse photon field. Making use of Einstein's equations (^, equation (^) can be recast 
as 



+ 



{dl - V^) A^ + 



/3e 



h<6h + ?,c + Ad)p-\{<oh + c)p 
o 2 



w^aT 



where Mp = (SttG) 



-1/2 



5. 



-i(66 - 3c - 2d)p + i(66 + 3c + 2rf)p ) (9^ A 
6 2 



T 



(102) 



2.4 X lO^^GeV is the reduced Planck mass. With H"^ = p/3M 



p = wp, and b = —5, c = 26 and d = —2, we finally arrive at 



1 + (26 + 33w;) J drA^ 



2/3e^(10 + 3^)] V^Af 



(103) 
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In section ^ we use this equation to study the time evolution of the photon field in infiation. 

In order to make a comparison with the photon coupling to a scalar field studied in 
Ref. [|15], first note that equation ( p.02| ) yields 



dr 



~ (l + ^/3e^[7-30e-(13-90eHj V^Af = 0, (104) 

where we made use of the Schwinger-DeWitt coefficients bg = (—5/16) + (15/8)^, Cg = 1/4 
and dg = —3/8, and Mg is the scalar field mass. This case is quite different from the fermionic 
case studied in section ^ in that for ^ < 1 /6 the photon field is damped in infiation. For a 
conformally coupled scalar field Eq. (|104|) reduces to 

d'rAr - (l + ^/?e^ (1 + ^)) = 0, (105) 

and no significant amplification occurs in inflation. On the other hand when ^ > 1/6 the 
photon coupled to a scalar field the amplitude grows in infiation. In this case the gauge field 
amplification can be obtained by effecting the following replacement in Eq. (H): 

^ = ^/3e (6e - 1) (1 - 3«;) e > ^ • (106) 

When both massive fermionic and scalar fields are present, one should add contributions 
from both fields. In de Sitter infiation [w = —1 and H = const) Eq. (|104|) simplifies to 

(l + f /5e^ (1 - 60) (5^ - V^) Ar = 0, (107) 
and thus conformal invariance is recovered, just like in the fermionic case (^. 

Appendix B 

We now outline computation of a volume averaged magnetic field strength for a given gauge 
field spectrum of the form 

Aj: = CoVh-'' cos{kT + ifio) , (108) 



where a = 6'^„ + 1/2 for the particular spectrum (]2l|) - (P2D, and we displayed explicitly the 



volume V dependence, which was for simplicity ignored in the Wronskian normalization ([T3|) . 
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Following Ref. |^ we define the volume averaged magnetic field correlated on a scale i as : 

B^, = {B,{i,x)B,{£,x)) - {B,{e,x)Bi{£,x)),^ , 



B,{i,x) = -— d'yB,{y), (109) 

4:7lt^ J\y~x\<e 

where the average (■) is taken over Fock space as well as the position x, and we subtracted 
the (divergent) vacuum contribution (■)vac- The definition ( |109| ) corresponds to a sharp ball 
window function w{x — y,i) = 1 for \x — y\ < 1, and zero otherwise. 

Ignoring for the moment the vacuum contribution, and making use of the Fourier trans- 
form 

B,{y, r)= j d'k'{27rfe"^'-yte,,ik^Ai{k', r) , Aiik', r) = ejApir) , (110) 

where ef denotes the unit transverse polarization 3- vector. To evaluate the incurring integrals 
it is convenient to introduce the following new variables : fi = y — x, f2 = z — x, K = 
{k' — k")/2 and k = {k' — k")/2, such that after some algebra we obtain 

/ Q \ 2 1 dk 
= (4^) 8V L ^^^'^^'^^ (2HCOS2H - sin2H)^ , (111) 

where we introduced the infrared and ultraviolet cut-offs, ko and A, respectively, which can 
be sent to infinity at the end of calculation. Keeping finite cutoffs is useful for discussion 
of the divergent limits, Ovac = 1/2 (vacuum spectrum) and Ogcinv. = 5/2 (scale-invariant 



spectrum). Equation (|111|) can be further simplified 



Bi = 5 Coi„, t„(xo,Xi)=/ .t^.d i [x) , (112) 



TT Jxn 



CO X ^ ' 

where xq = ^k^l, xi = 2Ai, ji = djo/dx = (cosx — sinx/x)/x denotes the spherical Bessel 
function, and we used the long time-average (cos^(A;r -|- fo))T = 1/2- When 1/2 < a < 5/2 
the integral La{xQ — > 0, xi — > 00) converges and can be evaluated exactly (c/. Eqs. 6.576.2 



and 9.122.1 in Ref. ||). The result is 



.,(0,00) = ^ V2 J V 2; ^ ^<a<^. (113) 



4 r(a)r(a + |) 
Putting everything together, we obtain 



5 9 X 22"-9 r f I - a) r fa - 1 5 

TT-^/^ T {a)T [a + 2 2 
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such that bi = y/S/ (S^/tt) for the thermal spectrum and 6| = S/Svr for the spectrum Bi oc i^^. 

For the vacuum spectrum (a = 1/2) the integral li diverges in the UV. Upon introducing 
an UV cut-off, the leading-log contribution to the integral equals ti = (1/2) ln2A£, and hence 
the spectrum reads 

Br = ^Co{ln2Aet^i-\ a = \. (115) 

The question is of course what is the natural value for the UV cut-off. At a first sight it may 
seem reasonable to choose the Hubble scale A ~ if, or the Planck scale A ~ Mp. We believe 
however that it is more natural to take A ~ e/2£, implying that measurements at some scale 
I are mostly sensitive to the field excitations at that scale. With this equation (|115D becomes 
cut-off independent : 

^ 87r\/2 " ' 2 ^ ' 

For the scale invariant spectrum a = 5/2, the integral La becomes infrared dominated, 
suggesting an infrared cut-off, ko. To leading-log accuracy ta(xo,oo) = — (l/9)lnxo {xq = 
2koi), so that 

i?;--- = ^Co (-ln2M)^ , a = ^. (117) 

In this case the natural cut-off is is the lowest momentum that gets amplified in inflation, 
implying that information about the beginning of inflation is present in the spectrum. This 
should not be taken very seriously however, because in practice it may be very difficult to 
distinguish a scale invariant spectrum from an almost scale invariant spectrum with a 5/2 
in which the dependence on the initial conditions in inflation drops out. For this reason, 
when plotting the scale-invariant spectrum we take /cq = e/2i, for which equation ( |117D 
becomes truly scale invariant, Bl"^ ™^- ~ Co/27r. 

Appendix C 

Here we consider two simple inflationary models: chaotic inflation and extended inflation. 
C.l Chaotic inflation with V{(f)) = mVV2 



For simplicity we assume that is the inflaton with the potential V{(f)) = m^0^/2 [p8 



where the cosmic microwave background radiation (CMBR) measurements constrain m ^ 
1.8 X 10^^ GeV. The considerations in this section can be easily generalized to any power-law 
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potential Vn{(f>) = {\n/n\){(f)/Mp)^, with A4 ~ 5 x 10 etc. The equation of motion for 
the inflaton reads 

d'^d) dd) r, , ^ 

-^+3ff-| + m=^ = 0. (118) 

where = p/3Mp, and p is the total energy density. Equation (|118| ) can be easily solved 



in the slow-roll approximation Q in which d'^(f)/dt'^ <^ 3Hd(j)/dt, d^V and {d(j)/dtY ^ 21^ 
such that 

i7^1^ = ^. (119) 
adt y/QMp 



The result is 



bo-J'^Mpmt, (120) 



where 0o can be related to the number of e- folds before the end of inflation, (pliN) 
4y/6Ml.N. Integrating ( p9|) once we then get for the scale factor 

a m(hnt m^t^ 



In 



ao VeMp 6 

Ht. (121) 



The time dependence of H becomes important only at late times in infiationj^ This then 
immediately implies 

(122) 



Ht 



and 

dd) dd) 2Ml 1 

-T = ~ _P_ 123 

dr dt (p T 

where ^^(r) ^ A^/QM^N . 
C.2 Extended inflation 

Extended inflation P0[ is realised by an exponential potential V{(f)) = Mpe~'^''^/^^^ , and it 



corresponds to the Jordan-Brans-Dicke theory of gravity rewritten in the Jordan frame. For 
an alternative formulation see [^. In its original version the model has limitations related 
to the graceful exit problem and the CMBR constraints. The spectral slope n of scalar 
•^For more formal treatment of the slow-roll conditions see for example 

^An alternative definition is the following time averaged expansion rate: H = Hq — w Hq — 



{ra^ I&Hq) ln(— iJoi"), where Hq = nup/^f&Mp. This definition gives the exact result for the scale factor (121) 
when expressed in terms of cosmic time t. 
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cosmological perturbations n — 1 = — A'^, and the fraction of perturbations in gravitational 
radiation r = 5A'^ are constrained to be much smaller than unity |^. Here we assume that 
the former problem is solved by matching at the end inflation onto a power law potential, 
while the latter can be resolved by an appropriate choice of the coupling constant A'. We 
discuss extended inflation primarily for its simplicity, and focus on the attractor solution |H2 



41| because it corresponds to more generic initial conditions than the slow-roll case. 
The infiaton equation of motion reads 

^ + 3H^ - A' Ae-AV/Mp ^ 0, (124) 
dt^ dt Mp ^ ^ 

where ^ 

4 g) (125) 

It is quite easy to show that equations (|124|-|125|) are solved by the following attractor solution: 



dt 



2Mp 

2Mp 
IT 



In 



Mp 



2 



such that 



and hence 



X'^ VA 



A'^ 



Ml 
t2 ' 



6 

A^ 



H 



1 da 
a dt 



Integrating this once gives 



(126) 

(127) 
(128) 
(129) 



a = {Mot f^' , 

where Mq < Mp is a mass scale that characterises the beginning of inflation, to ~ ^/Mq. 
This immediately implies that A' < \/2 is required to get a superluminal expansion. On the 
other hand to get a nearly scale-invariant spectrum of perturbations. A' ^ 1. Since inflation 
is assumed to terminate at ~ 0, in order to get more than 60 e- folds (required to solve 
the standard cosmological problems), the potential Vq should be sufficiently small, such that 
at the beginning of inflation we have V(0o) Mp. Now making use of adr = dt we can 
reexpress ( |129D as 

-llMor) , (130) 



A' 



36 



so that 



Finally, equation 



dr dt 



X' 



2Mf 



(131) 



can be expressed as a function of conformal time as follows: 



2A' , 
^Mp In 

2 - 



p - 1 ) MoT 



In 



/_6_ _ \ M|,Mg 



(132) 



Equations (|123| ) and ( |131| ) are used in section ^ for computation of magnetic field amplifi- 
cation in inflation. 
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